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1. Introduction 

We investigate the local and global bifurcations of bulk fluid flow to determine its 
dynamical state. These flows occur during crystal growth of mercurous chloride by 
the physical vapor transport (PVT) process. The dynamical state of the flow field is 
analyzed for a parameter space which corresponds to ground based experimental condi- 
tions. In particular, two sets of experiments are considered which corresponds to: (a) 
Laser Doppler Velocimetry (LDV) measurements of the flow field by Kim et al. [1], 
and (b) Correlation of crystal quality with growth conditions of mercurous chloride, 
Singh et al. [2]. Crystal quality is correlated to the scattering patterns of light through 
the crystals. In view of these experiments, we analyze the physics of PVT inside a rec- 
tangular enclosure through a mathematical model to provide insight into various trends. 
A parametric setting is selected which will lead to a meaningful microgravity experi- 
ment 

PVT is advantageous for a material such as mercurous chloride due to its decom- 
position near its melting point [3], and its high vapor pressure at typical growth condi- 
tions (T = 345 °C ,P A =289 torr). Various vapor transport processes and their convective 
characteristics are discussed by Rosenberger [4]. The wide infra-red transmission range 
(.36 to 20 \im) of mercurous chloride makes it attractive for applications in acousto- 
optic optoelectronic devices such as acousto-optic tunable filters, Bragg cells, and spec- 
trum analyzers. PVT involves sublimation from a heated source material, transport 
through a rarified inert gas, and condensation to a lower temperature heat sink inside 
an enclosure. The dynamical state of the flow field during PVT is affected by the 
external conditions through the thermal settings, gaseous components, ampoule 
geometry, and orientation of the enclosure with respect to gravity. 



There have been numerous works which investigate PVT under various conditions 
and ampoule geometry. Greenweil, Markham and Rosenberger [5] show that, even in 
the absence of gravity, non-uniform concentration gradients occur near the vapor cry- 
stal interface for cylindrical ampoules oriented horizontally. For vertically oriented 
cylinders, Markham Greenweil and Rosenberger [6] show that convection is stabilized 
for solutally/thermally dominated flow when transport is against/assisted by gravity; 
otherwise convective cells occur in both systems beyond a certain parametric thres- 
hold. The same trends are shown to occur for a square enclosure oriented vertically 
[7,8], which also agrees with the results of Markham and Rosenberger [9]. Effects of 
surface reactions, and irregular thermal profiles are considered by Zappoli [10] and 
Extremet et al. [11] for horizontal rectangular enclosures. Zhou and Zebib [12] show 
the parametric range for which flow fields bifurcate from steady to oscillatory for dou- 
ble diffusive flows inside horizontal enclosures in PVT. This model considers PVT 
inside a vertical enclosure for which the dynamical state of the flow field is beyond the 
critical state for a range of parameters. 

The purpose of this research is to examine computationally the local and global 
bifurcations of flow fields in PVT for practical ground-based crystal growth conditions 
in order to predict their dynamical states, and to delineate growth conditions for a 
microgravity experiment. This will allow testing of the hypothesis that for a system far 
from thermodynamic equilibrium, a change in the gravitational level can shift its 
dynamical state. The results are compared to LDV measurements, and other experi- 
ments to correlate crystal quality with practical growth conditions. The local bifurca- 
tions show that for a fixed aspect ratio, as the Rayleigh number increases the flow 
bifurcates from steady to oscillatory for practical growth conditions. The steady regime 
consists of both diffusive-advective and diffusive-convective flows. The global bifurca- 
tion of the flow shows that in the neighborhood of the critical Rayleigh number the 
flow field bifurcates from uniform flow to an asymmetric cell. This corresponds to the 
transition from diffusive-advective to diffusive-convective flow. Subsequent global 
bifurcations as the Rayleigh number increases show a bifurcation to a symmetric dou- 
ble cell structure. Depending on the range of the Rayleigh number, the symmetric 
structure bifurcates to a number of asymmetric modes which leads to a traveling wave. 
The traveling wave either decays or persists. In the nonlinear regime, the flow switches 
from steady to aperiodic. Computations of asymptotic states of the flow show that 
these oscillations do not decay and proceed during the entire crystal growth. We show 
the existence of a traveling wave for this parametric range. Traveling wave state has 
been shown to occur in double-diffusive Rayleigh-Befiard systems, Yahata [13]. The 
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traveling wave state is correlated with a decrease in crystal quality [2] as measured 
through light scattering for the given parametric value. We show, for the condition 
which corresponds to the traveling wave state, that a typical microgravity environment 
of icr 3 g 0 stabilizes the flow field and results in a diffusive-advective flow. 

This paper is organized as follows. We review briefly the physics of PVT and 
then compare our computational results to recent LDV experiments. The dynamical 
state of the flow field is correlated with crystal quality for a range of practical crystal 
growth conditions. From the analysis of these results a critical microgravity experiment 
is proposed. 

2. Mathematical Formulation 

PVT occurs inside a vertical enclosure for source temperature T s greater than the 
sink crystal temperature T c , Figure 1. For AT > 0, the source material A sublimes 
through an inert component B, and condenses at the sink. The inert component B is 
assumed to be rejected at the source/sink interfaces. The phase boundaries at the 
source and sink are at thermodynamic equilibrium, which fixes the mass fractions, co s 
and coc . For a system which deviates from equilibrium, AT > 0, motion of the binary 
components inside the enclosure obeys the conservation of mass, energy, balance of 
momentum, and species continuity. Restricting the binary system to thermal convec- 
tion, M a -M b , the governing equations for the average (barycentric) velocity are res- 
tricted to the Boussinesq equations as follows: 


V-V* - 0 

(1) 

. - - Vp' +jlv 2 v* + pf 

(2) 

El - a v 2 r* 

(3) 
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with the following initial and boundary conditions, 
Initial conditions: 


t - 0 0 < x* < L , 0 < y' < H 

T*(x*,/,0; )-r s * + (y*IH) (T c -T s ) (5a) 

©a(** j’,0) - (o^s + (y* IH) ((o AC - ©as) (5b) 

Boundary conditions: 

No slip at the left and right walls, 

t > 0 v* ft - 0 , x* - 0 , X* - H , V y (6) 

Prescribed temperature and fixed mass fraction at the source and crystal boundary, 

T'(x* ,0,t') = T s , ©a(x’, 0, r *) -©as , y*-0 , V x (7a) 

T'(x\H,t*)-T c , a>A*(x*,//,r*)-o>Ac , / - H , V x (7b) 


Insulated and impermeable boundaries at left and right walls, 

Vr* ft — 0 , Vo>a • ft — 0 on T at x* — 0 , x* — L , V y (8) 

Mass flux condition at the boundary of the source and crystal, 

V* ft-- , D *f . (V©A - ft) , at y" - 0 , y* -H , V x (9) 

1 - © a (* , 0 , t ) 

with vanishing tangential velocity at the interface, 

V* ?-0 on r at y* -0, y' -H , V x (10) 

In the above equations, the normal to the internal boundary is denoted by ft, t its 
tangent, and r its boundary. The starred (*) variables are dimensional, V ■ for all , and 
the overbar represents average values. The mass flux condition which consists of 
Fickian diffusion and bulk advective flow is relative to the fixed crystal boundary. This 
condition was derived by Westphal and Rosenberger [14]. The motion of the source 

and crystal is assumed negligible, since the velocity of the bulk fluid is much greater 

than the crystal growth rate. We selected to use insulating boundaries to emphasize the 
effects of convection. Studies comparing effects of conducting versus ins ulating boun- 
daries for confined fluids inside enclosures [15] show insulating boundaries are less 
convectively stabilizing Ra c ins < Ra c cond . Thus the critical value for insulating boun- 
daries represents a conservative estimate. 
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PVT occurs at constant total pressure, 

Pt “Pa + Pb 


The partial pressure for mercurous-chloride is estimated using the correlation 
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with a =29.75 and b = 11767.1 yield P A in Pascal units [16]. Since the presence of 
imp urity gases is inevitable, purification experiments [17] indicate the likelihood of the 
presence of CO. We assume P B - 5 torr for the estimation of the thermophysical pro- 
perties from the Chapman Enskog relationships [18]. The presence of the impurity 
component is used as a perturbation to the thermophysical properties of mercurous 
chloride, this allows for a more realistic estimate of the parametric values. The gas 
mixture is assumed to be ideal with 


P-1 IT 


and 


(o A = 


1 

1 + WTJPa - 1 ) m b im a ' 


(13) 


using the mass and mole fraction relations. 

To determine the dynamical state of the flow field requires the solution 
? (x ,;y,z,r,n) as a function of the 3 space coordinates (x,y,z), time it), and a vector 
space n of the independent parameters. Even with the present computational capabil- 
ity, the solution of the entire space-time continuum (4-dimensional) is a formidable 
task. However, for a subset of the parameter space n, it is possible to reduce this task 
and find 3-dimensional steady state solutions V (x ,y ^ ;I1). This approach will yield ade- 
quate solutions for parametric values up to the Hopf bifurcation, which is a transition 
from the steady to periodic state. Even though this approach is desirable for many cir- 
cumstances, limitations exist. Estimation of the parameter space n indicates that the 
critical parameters are outside the range for which steady-state solutions exist. Alterna- 
tively the 3-dimensional solutions of V* (x,y,f;II) can be obtained. This is commonly 
referred to as 2-dimensional unsteady, even though time is a valid dimension. To study 
the dynamical state of the flow field, we selected to find solutions for V (x ,v ,t ;II). The 
advantage to this approach is that the state of flow field can be predicted for a much 
wider range of the parameters in n. However, Bontoux et al. [19] has shown that the 
critical parameters for 3-D steady state models are less than those predicted by 2-D 
models. Thus the critical values for our model represent, at best, conservative estimates 
for the solution of the entire space-time continuum. Even with the inherent limitations 
of our model, it serves as a tool to study the fluid physics of PVT beyond the steady- 
state solutions. 
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2.1 Scaling 

The restriction of the above equations to two spatial dimensions reduces the 
momentum equations to a single equation for the vorticity field, and allows for com- 
pact representation to vorticity-stream function relationship for computation. The equa- 
tions are scaled using the characteristic time, and velocity (t/ c ) as L 2 la, and aJL respec- 
tively. Independent characteristic lengths L and H are used for the horizontal and verti- 
cal directions. The temperature and mass fraction are normalized as 7" - (T* - Tq)IAI, 
o) A = (©; - coJJ c )/A<b, in which AT and Aco equal (Ts - Tc ) and - a> AC ), respectively. 
Let, 


« d'P* > _d x ¥* dv* _ du* 
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(14) 


The above scales yield the following dimensionless equations: 
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(16) 

(17) 

(18) 


The mass flux conditions at the boundaries become: 

Ar 1 


v(x,0,t) - - 
v(x,l,t) - - 


Te (Cv — 1) dy 

Ar 1 da>AC 
Te Cv dy 


(19a) 

(19b) 


The dimensionless parameters are: 


-4 


Ra 


PAT ng 0 L 3 


Pr - — , Le 
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va 


a 

&ab 


Cv 


Dab 

(i - ©; c ) 
Sa> 


The scaling gives rise to six independent parameters. The vector space for the 
parameters may be expressed as 


II - II (Ra , Ar , Pe , Pr , Le , Cv) 


( 20 ) 
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These parameters are respectively, the Rayleigh number (Ra), the aspect ratio (Ar), the 
Peclet (Pe), Prandtl (Pr), Lewis (Le) numbers, and the concentration parameter (Cv). 
Note the use of multiple length scales, L & H for the horizontal and vertical direction, 
result in the explicit appearance of Ar in the governing equations. In order to show 
the explicit dependence of the mass flux condition on the Peclet number, we used a 
characteristic velocity U 0 which is independent of the velocity U c due to convection. 
This condition occurs because the base-state flow field is independent of buoyancy 
induced convection. Even if g - 0, there always exist a base diffusive-advective flow 
for AT > 0 which is dependent on the thermodynamics. This dependence has been 
shown [5] through the relationship of Pe on the partial pressure as 

Pe - ln(P B (H)IP b (0)) . (21) 

This equation can also be expressed as a function of the mole fraction, via, Xi = Pi IPt • 
Thus the value of Pe is fixed by the thermodynamics of the process. 

The above equations are solved using finite difference techniques [7]. For the 
aspect ratio range considered, we used a 40x90 grid size. In addition a third-order 
Adams-Bashforth technique was used for time discretization. The choice of grid size 
and time step was adequate to resolve the details of the flow field. 

3. Results and Discussion 

3.1 Parametric Space Variation 

The parametric range for which we investigated the dynamical state of the flow 
field is shown in Table 1. Cases 1-3 correspond to experimental conditions for the 
recent LDV measurements of Kim et al. [1] of the radial velocity profile at a fixed 
axial location. Whereas, cases 4-9 correspond to crystal growth experiments of Singh 
et al. [2], which correlated crystal quality with crystal growth conditions by varying 
the source and sink temperature. Case 10 corresponds to case 8 in a microgravity 
environment with a nominal gravitational acceleration of l(r 3 g 0 . These experiments 
were conducted in cylindrical ampoules; we use dynamic similarity [20] for application 
in rectangular coordinates in which the width of our enclosure is set equal to the diam- 
eter (L = D ) for our parametric investigation. The aspect ratio for the cases considered 
is on the order of .12, the Rayleigh number range from 2.08X10 2 to 9.93x10 s . 

3.2 Local Bifurcation 

Our parameter space n dictates a functional space of six parameters. However, 
for the parametric range of interest, most of these parameters vary over a narrow 
range and can be considered as passive. This reduces the parametric space to a 



codimension-one bifurcation n - n(fta). The local bifurcation of the flow field to del- 
ineate between steady (stationary state) and oscillatory convection (periodic state) is 
shown in Figure 2 as a function of the Rayleigh number. In the steady convective 
region the flow field may be either diffusive-advective or diffusive-convective. The 
diffusive-advective flow is one-dimensional and is also known as the Stefan wind. It 
corresponds to the base flow in PVT by which global bifurcations lead to various 
modes as the Rayleigh number increases. If the flow field becomes unstable as the Ra 
increases, the first threshold of instability leads to the diffusive-convective regime, 
whereby a number of modes can occur. Transient oscillations can also occur in this 
region, however, long time asymptotic computations show that these oscillations 
decay. The second bifurcation connect the stationary state to the periodic state and 
occurs in the neighborhood of 2.59x10* < Ra £2.08xl0 5 . In this region the flow field 
oscillations do not decay. Time asymptotic computations indicate that these oscillations 
proceed for the entire growth of the crystal. This periodic state causes unsteady heat 
and mass transport which affect the homogeneity in refractive index of mercurous 
chloride crystals, and hence results in a degradation of crystal quality. 

3.3 Global Bifurcations 

For the Rayleigh number range in Figure 2 beyond the diffusive-advective region, 
there exists global bifurcations of the flow field which lead to various flow field plan- 
forms. This is illustrated for cases 1-2 in Figures 3-4, case 3 is similar to case 2. The 
evolution of the flow field in Figure 3, case 1 Ra - 1.01x10*, shows that the flow 
becomes unstable against the base diffusive-advective flow (t - 0.34). This instability 
leads to an asymmetric cell (t - 0.66) which grows until it occupies the entire cavity 
( t = 1.04). The asymmetric unicell is further bifurcated into two asymmetric cells 
(t - 151). The long-time asymptotic solution shows further bifurcation to an asym- 
metric cell which is pushed toward the crystal. This asymmetric structure is sustained 
for cases 2 & 3, however the flow becomes diffusive-advective for case 1. 

The difference between the flow field in vertical PVT and typical Rayleigh- 
Benard (RB) convection in cavities is the mass flux through the source and sink boun- 
daries quantified by the Peclet number. In addition in RB systems the first threshold of 
instability is against a quiescent state in comparison to the diffusive-advective flow in 
our case. The mass flux through the boundaries lifts the threshold of instability, thus 
Figure 2 for a RB system would be shifted to the left. The similarity between the two 
systems is exemplified in the experiments of Olson and Rosenberger [21] to study ins- 
tability of flow fields in bottom heated vertical cylinders. Their experiment showed that 
the first threshold of instability leads to an asymmetric roll; the number of rolls 
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increase with Rayleigh number. This convection planform is in agreement with our 
results for rectangular geometry. In contrast, the 2-D axisymmetric model of Kim et al. 
[22] for PVT in cylindrical geometry predict symmetric cells. This is due to the 
axisymmetric constraint imposed in the 2-D model. The ability of a rectangular model 
to predict the convective planform in cylindrical geometry underscores the importance 
of dynamical similarity to investigate trends. 

An increase in the Rayleigh number (Ra = 2.77X10 4 ) shows that the global bifurca- 
tion from diffusive-advective flow leads to symmetric cells (t « 0.19). This is in agree- 
ment with the experimental findings of Olson and Rosenberger [21] for the second 
bifurcation. The symmetry is not sustained and bifurcates to a superposed flow consist- 
ing of an asymmetric cell and uniform flow near the left wall (t - 0.31). The flow field 
is subsequently transformed from a core driven flow to a boundary layer driven flow 
indicated by the temperature and concentration contours (t - 0.40); at this point the flow 
field is dominated by the asymmetric cell. This asymmetric cell is subsequently bifur- 
cated to a number of asymmetric cells (/ - 0.44). Subsequent bifurcations t 2: 0.52 show 
the existence of an extended traveling wave. The long time asymptotic solution shows 
that the traveling wave decays, and the flow becomes steady with a single cell asym- 
metric planform. The time history of the flow field at a fixed point in the flow also 
shows the decay of the oscillations. 

3.4 Comparison with LDV Experiments 

Case 2 corresponds to the experimental condition for which the velocity profile 
was measured by LDV. In spite of the limitations of our model as we had pointed out, 
comparison with the experimental data of Kim et al. [1] show reasonable agreement. 
Figure 5. The magnitude of the flow is within the experimental measurements. How- 
ever, the model shows a steeper velocity gradient near the wall. The asymmetry of the 
velocity profile is correctly predicted by our model. Case 3 shows similar agreement, 
however, the direction of the flow field is opposite to that predicted by the experi- 
ments. This is reminescent of Rayleigh-Benard cells inside a cavity for which equal 
probabilities exist for the sense of rotation of the cells. Case 1 also shows reasonable 
agreement for short time, however, long time asymptotic computation shows that the 
flow field gains stability. In view of the recent work by Rosenberger et al. [23], it 
must be pointed out that even though this 2-D unsteady model predicts reasonable 
trends in the core of the flow field, it may not be adequate for the prediction of crystal 
growth rates. 

3.5 Correlation of Crystal Quality with Dynamical State of the Flow Field 
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To determine the effects of crystal growth experimental conditions, hence convec- 
tion, on crystal quality, a number of experiments were carried out by Singh et al. [2], 
The temperature difference between source and sink were varied in order to effectively 
vary the transport conditions which are characterized by the Rayleigh number and the 
re maining parameters, cases 4-9 in Table 1. The crystal quality was quantified by light 
scattering through the crystal, birefringence interferograms, and x-ray rocking curves. 
All three techniques indicated a deterioration in crystal quality with increasing Ray- 
leigh number, cases 5,6,8. An example of the light scattering results as a function of 
the Rayleigh number is shown in Figure 6. The spread in the sidelobes for the largest 
Rayleigh number indicates an increase in inhomogeneity which causes variation in the 
refractive index. The inhomogeneity is attributed to the convective state during crystal 
growth, since all other parameters were kept fixed. 

The convective characteristics as a function of the Rayleigh number, for the cry- 
stals in which the quality was quantified, is shown in Figure 7. For the lowest Ray- 
leigh number ( Ra - 6.48X10 2 ), the flow is diffusive-advective and represents an optimal 
condition for crystal growth. As shown in Figure 6, the best quality crystal was 
obtained for those conditions. For increasing Rayleigh number, a sequence of bifurca- 
tions occur which results in various asymmetric modes. The global bifurcations for 
cases 6 & 7 are similar to cases 1 & 2, thus the time sequences are not shown. For 
Ra « 1.06X10 4 , the first bifurcation from the diffusive-advective flow yields an asym- 
metric cell which is further bifurcated to two cells as shown. In contrast, for 
Ra = 2.59x10* the first bifurcation yields two symmetric cells, similarly to case 2 tran- 
sient oscillatory convection which results in a traveling wave occurs. Asymptotic solu- 
tions show that the flow field for both cases 6 & 7 become steady and yields an asym- 
metric cell near the crystal. However, for for Ra s 2.08x10 s , the oscillations do not 
decay in time. The global bifurcation sequences illustrating the traveling wave state for 
this case are shown in Figure 8. The flow bifurcates from diffusive-advective to a sym- 
metric mode and subsequently to multiple asymmetric cells (/ = 0.059). As time 
increases an extended traveling wave propagates alternatively from source to sink 
t - 0.081 -» 0.177. Defect structures are also formed during the wave propagation. Long 
time asymptotic solutions show a transition from extended to localized traveling wave 
near the crystal interface. The traveling wave state is further illustrated in Figure 9 by 
the stream surface plots. Propagation of the traveling wave near the crystal and its 
reflection as it travels toward the source is shown. The dynamics for case 9 is similar 
to case 8. 
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4. Dynamical State of the Flow Field 

The dynamical state of the flow field for the parametric range in Table 1 shows 
that the state of the flow field range from steady to aperiodic. In some cases the oscil- 
lations decay, in order insure that the dynamical state of the flow field is correctly 
predicted, we carried out asymptotic computations as pointed out by Goodrich et al. 
[24] for the representative cases. Large storage requirements are necessary for these 
cases on the order of 1 gigabyte of data. The dynamical state is predicted from the 
time history of the flow, phase space trajectory and power spectrum. We use the 
Cooley-Tukey algorithm [25] with fast fourier transforms to evaluate the power spec- 
trum. The time step is small, such that the Nyquist frequency criterion is satisfied, thus 
eliminating aliasing. The power spectrum and its smoothing is calculated respectively 
from 


P u (f) = {*-} e ' i2 ' f ‘ dt 2 

(22) 

Py(f)~W S (f)* P v (f) ■ 

(23) 


W s (f) is a Hanning spectral window used in the convolution relationship with the 
estimated power spectrum P v (f) to obtain the smooth spectral estimate. The time 
interval is denoted by R. 

The dynamical state of the flow field is illustrated in Figure 10. The time history 
for case 5, Ra = 3.69X10 3 , which corresponds to diffusive-advective flow shows an 
asymptotic approach to steady state. In contrast, in the diffusive-convective region, 
case 7, the flow field undergoes an aperiodic transient before settling to a steady state. 
An example of the bifurcation from steady to an aperiodic state is shown for 
Ra = 2.08x10 s , case 8. Long time asymptotic computations show that the fluctuations 
never decay, they proceed during the entire crystal growth event The broadband 
features of the power spectrum and the unpredictability of the system illustrated by the 
thick limit set in phase space indicate a chaotic state. Thus, it is not surprising that 
lower quality crystals result for this flow state. However, conducting case 8, in a 
microgravity environment of l0' 3 g o , case 10, shows that the flow field becomes 
diffusive-advective similar to case 5 in Figure 7. An alternative method to stabilize 
the flow field is the use of narrow aspect ratio cavities, however, a large decrease in 
size of the crystals may not be suitable for device applications. Thus, the parametric 
range for high quality crystal growth on earth is limited. In particular, stable flow 
fields might not exist as the aspect ratio increases. In addition, as we had pointed out, 
in 3-D the critical values are lower, thus the bifurcations would be shifted to the left in 
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Figure 2. To obtain maximum benefit from a microgravity environment, a parametric 
range for which the dynamical state of the flow field is aperiodic should be selected as 
in case 8 for a microgravity experiment. This condition would shift the dynamical state 
of the flow field from chaotic to steady. Thereby, the hypothesis that for a system far 
from thermodynamic equilibrium, a change in the gravitational level can shift its 
dynamical state, can be properly tested. 

5. Conclusions 

We examined the local and global bifurcations of the flow field in PVT for a 
parametric range for which experiments, to measure the velocity profile, and to quan- 
tify the effect of these parameters on crystal quality, have been performed. The phy- 
sics of PVT as predicted by our model shows reasonable agreement with the velocity 
profile measurements. Correlation of crystal quality with the fluid dynamic state of the 
flow field shows that a chaotic dynamical state is correlated with lower crystal quality. 
In this chaotic state, a traveling wave propagates in the flow field and provides the 
mechanism for unsteady transport However, a microgravity environment, with a g- 
level of ltr 3 g 0 , effectively stabilizes the flow field and provides an optimal condition 
for high quality crystal growth. In essence, the hypothesis that for a system far from 
thermodynamic equilibrium, a change in the gravitational level can shift its dynamical 
state, has been shown theoretically. The dynamical state of the flow shifted from 
chaotic to steady. The results from a microgravity experiment would prove the 
hypothesis. 
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7. Nomenclature 

Ar aspect ratio (L/H) 

Dab molecular mass diffusion coefficient 
g 0 acceleration of gravity on earth (9.8 ml sec 1 ) 

H cavity height 

Hm axial location of radial velocity profile measurement 
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f unit vector in vertical direction 

L cavity width 

Le Lewis number 

M molecular weight 

n ratio denoting reduction in acceleration of gravity 

p hydrodynamic pressure 

P partial pressure 

p v power spectrum of variable v 

Pe Peclect number based on mass diffusion 

Pr Prandtl number 

R time duration of data 

Ra T thermal Rayleigh number 

t dimensionless time 

T dimensionless temperature 

T, characteristic time (L 2 /a) 

u, v dimensionless velocity in x and y directions (u* IU C ,u*IU c ) 
U c characteristic velocity based on convection (a/L) 

U 0 characteristic velocity based on the advective-diffusive flux 

x, y dimensionless horizontal and vertical directions (x*IL ,y*IH ) 
W s spectral window 

V for all 

Greek characters 

r boundary of cavity 

n parameter space 

a thermal diffusivity 

p thermal expansivity 

A difference 
X mole fraction 

|i dynamic viscosity 

v kinematic viscosity 

to mass fraction 

4* dimensionless stream function (¥’ 1U C L) 

V gradient in x and y directions 



p density 

V 2 Laplacian in two directions 

£ dimensionless vorticity LIU c ) 


Subscripts 

A, B,i denote components 
c crystal 
s source 
T total 

Superscripts 

* dimensional quantity 
- average quantity 
c critical 
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Table 1 - Case Studies Corresponding to Velocimetry and light Scattering Experiments 


Case 

T S (C) 

T c (C) 

AT 0 

L(cm) 

H (cm) 

Ar 

Ra 

Pr 

Le 

Pe 

Cv 

i 

310 

290 

20 

22 

17.7 



.899 

.404 

2.49 

1.091 



290 

30 

22 

17.7 

.124 

2.77 X 10 4 

.894 

364 

3.03 

1.051 



290 

20 

2.3 

18.0 

.127 

1.15 X 10 4 

.899 

.404 

2.49 

1.091 

H 

327 

323 

a 

1.1 

9.0 

.122 

6.48 X 10 2 

.911 

.392 

1.74 

1.213 

5 

346 

338 

8 

1.1 






2.79 

1.065 

6 

360 

349 

11 

u 

9.0 





3.42 

1.034 

■ 

373 

359 

14 

1.1 

9.0 

.122 

2 J 9 X 10 4 

.907 

.392 

3.95 

1.019 

n 


390 

20 

1.1 

9.0 

.122 


.905 


5.13 

1.006 

n 


409 

29 

1.1 

9.0 

.122 




6.02 

1.002 

+10 

410 

390 

20 

1.1 





.390 

5.13 

1.006 


+ Microgravity condition 10“ 3 g o corresponding to case 8 
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T s > T c 

w s > w c 

A - crystal component 
B - inert component 
(T s = 345 °C, P A = 289 torr) 


Figure 1 .—Description of physical vapor transport 
process at typical operating conditions. 
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Figure 6. — Light scattering patterns through crystals as 
a function of Rayleigh number after Singh et al. [2]. 
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